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§ ; Abstract 

'The transverse electron scattering response function of ^He is studied in the quasi-elastic peak region for 
^^^nomentum transfers between 500 and 700 MeV/c. A conventional description of the process leads to results at a 
^Substantial variation with experiment. To improve the results, the present calculation is done in a reference frame 
Tthe ANB or Active Nucleon Breit frame) which diminishes the influence of relativistic effects on nuclear states. 

o 

^^i'he laboratory frame response function is then obtained via a kinematics transformation. In addition, a one- 
^^pdy nuclear current operator is employed that includes all leading order relativistic corrections. Multipoles of 
^^his operator are listed. It is shown that the use of the ANB frame leads to a sizable shift of the quasi-elastic peak 
^o lower energy and, contrary to the relativistic current, also to an increase of the peak height. The additionally 
considered meson exchange current contribution is quite small in the peak region. In comparison with experiment 
one finds an excellent agreement of the peak positions. The peak height agrees well with experiment for the lowest 
considered momentum transfer (500 MeV/c), but tends to be too high for higher momentum transfer (10% at 700 
MeV/c). 

PACS numbers: 25.30.Fj, 21.45.-v, 21.30.-x 



1 



I. INTRODUCTION 



we studied the longitudinal electron scattering response function of trinucleons. We, as well 

m 

observed that for increasing momentum momentum transfer g, in particular for q > 500 



In Ref . 

as others, 

MeV/c, the non-relativistic theoretical results increasingly deviate from experiment. A similar problem 
arises in the case of the transverse response of trinucleons {2, 3, 4]. These problems appear to be related 
in part to a deficiency of the non-relativistic nuclear dynamics at such q values. In Ref. [jSi] methods were 
proposed which would allow the extension of such non-relativistic calculations to higher q. These methods 
proved to be efficient in the case of the longitudinal response. 

In the present work, with the help of such a method we analyse the transverse response function of ^He 
in the quasi-elastic peak region. Another improvement on the non-relativistic description in the present 
work results from our taking into account all the leading order relativistic corrections to the one-body 
electromagnetic current operator. Such corrections have been employed in the deuteron case js] and they 
were included {2, 2] when calculating magnetic form factors for elastic electron scattering on trinucleons. 
However they have not been previously taken into account for the A=3 transverse responses. Here we 
account for these corrections via considering the current operator that contains all the correction terms 



of the M order. We calculate this operator proceeding from corresponding matrix elements jsl of the 
current. 

Our preceding study of the transverse response of He [4] was done in the framework of a non-relativistic 
description with inclusion of full final state interaction via the Lorentz integral transform method [l, 8]. 
We used the BonnA NN potential [o^ plus conventional NNN forces as a nuclear dynamics input. Use of 
the BonnA potential gave a "unique" prescription for meson exchange contributions to the electromagnetic 
current of the nucleus. It is of interest, however, to use for our present study more modern NN interactions 
such as the AV18 potential [10]. Our results for the transverse response functions of the trinucleons using 



the AVIS NN plus UIX NNN 



ll| potentials have recently appeared 1^ for the threshold region. In that 



13| meson exchange currents 



paper we describe our procedure of using the Arenhovel-Schwamb technique 
of the AVI 8 potential. In the present work we extend our considerations to the quasi-elastic region and 
consider various intermediate momentum transfers. 
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II. FORMULATION 



In the one photon exchange approximation the cross section for the process of inclusive electron scat- 
tering on a nucleus is given by 



_fa_ 
dfl duj 



(^Mott 



-RL(giab,t^lab) + 7^ +tan^- 1 i?T(giab,C^L 



abJ 



^fab \^ 2giab 2 

where and Rt are the longitudinal and transverse response functions respectively, cuiab is the electron 
energy loss, giab is the magnitude of the electron momentum transfer, 9 is the electron scattering angle, 
and = g2^^ _ ^2^^. 

In the present work we study the transverse response function. It may be written down as 



Rriqi^h, wiab) = ^ dfi^t)if ■ {^t)fi5{Ej - Ei - o^iab) 



(2) 



Here the subscripts i and / label, respectively, an initial state and final states, including their total 
momenta P, and Pj. One may write df = dPjdf. Eq. ([2]) contains df only. The notation Ei, Ej refers 
to total initial and final-state energies. (In jj] the notation Eij was used for internal energies.) The 
quantities ( J^) are on-shell matrix elements of the transverse component of the nuclear current operator 
J(q,^), 

(Ji)/,5(P/-P,-q) = (vl/^-|J<(q,^)l^.), (3) 

taken at q = qiab, ^ = i^iab, Pj = 0. The states entering here are eigenstates of the total Hamiltonian 
with eigenenergies Ef and Ei. They are normalized as 

(*/|*/,) =5(/-f), {^^\^e) = S{P,-P,). (4) 

The above relationships refer to the laboratory reference frame. It is useful also to consider a response- 
type quantity defined by the same relationships referring to another reference frame. We shall denote 
the corresponding quantities g^, etc. In particular, the states and then will be eigenstates of the 
total Hamiltonian in the reference frame considered. For the class of reference frames moving with respect 
to the laboratory frame along the q direction the following relationship is valid: 

-Rr(9iab, i^iab) = T-^-Rrlfc; '^fr)- (5) 

IVIjn 

Here Mj- is the mass of the target. 

Relativistic effects are present in Eq. ([3]) both in the states and in the nuclear current operator. 

To account for the former effects we proceed as in the longitudinal case [Sj and introduce the active nucleon 



Breit (ANB) frame. In the ANB frame, the nucleus has the momentum — y4qANB/2 in the initial state, 
Qanb being the momentum transfer from the electron to the nucleus in this reference frame. At high q 
values, nucleon momenta in the initial state have the values of about — qANB/2 in this reference frame. 
In the final state in quasi- free kinematics the active nucleon has a momentum about qANB/2 while the 
momentum of each of the other nucleons remains at about —qANB/2. Thus, typical initial and final state 
nucleon momenta are restricted to magnitudes of about (3'anb/2 — q/2 in the ANB reference frame while, 
say, in the laboratory frame nucleon momenta up to q are present. Furthermore, it also follows from 
the above that the energy transfer ciJanb in the ANB reference frame is zero at the quasi-elastic peak, 
and this applies both to the relativistic and the non-relativistic case. Therefore, even when one treats 
the nucleus non-relativistically the peak remains at the same position as in the relativistic case. This 
contrasts with a description of the process in the laboratory reference frame where positions of the peak in 
the relativistic and the non-relativistic cases would differ considerably. Hence non-relativistic calculations 
in the quasi-elastic region should be done in the ANB frame to minimize errors due to relativistic effects. 
The laboratory response function sought for is obtained subsequently with the help of Eq. ([5]) with "ANB" 
being substituted for "fr". 

We perform the corresponding non-relativistic calculation in the ANB reference frame. One defines 
the internal current operator J obtained by taking a matrix element in the center of mass subspace of the 
total current operator: 

J5(P/ - P. - q) = (P/|J(q,o;)|P,). (6) 
At Pj = Pjq, q being g^^q, this operator may be written as J(q, uo, Pi). One may then rewrite Eq. ([2]) as 

Rt^^ {q ANB, ^anb) = ^ df{i!i\3l\i!f) ■ {ipf\Jt\ipi)S (e/ - e(gANB, luanb)) , (7) 

where the transverse component Jt of J(q, co',Pj) is used and the values q = qANB, ^ = i^anb, 
Pi = — ylgANB/2 are set. Here ipi and ipf are the non-relativistic internal states. They are independent of 
the center of mass momenta. The energy e/ is the internal energy in the final state, and 

(^ANB)^_(pANB)^ iA-l)qU 

e(,gANB, i^anbJ = + cjanb H ^ttt = Cj -h cuanb H ^ttt , 

Zlvij' ZIVl'jp 

where is the internal energy in the initial state. One also has 

gANB = 7(?lab - ^C^lab), C^ANB = 7(^lab " /^^lab) , 7 = (1 - 



2{Mt/A) [ 2{Mt/A) 



(One gets wanb = when substituting cjiab = [{Mt/AY + q^^^ ' — Mt/ A in the expression for cuanb- 
This is in agreement with the said above.) 



III. THE NUCLEAR CURRENT OPERATOR AND ITS MULTIPOLE DECOMPOSITION 



We employ the transition current operator that is a sum of one-body and two-body currents. In [4| 
the non-relativistic expression for the one-body current was used. For the present apphcations we have 
calculated relativistic corrections to the one-body current operator. To do this we proceeded from the 
expressions for matrix elements of the one-body current of the form (p/|J|pj) listed in Ref. p]. The 
operator so obtained reproduces these expressions. 

This operator denoted J'-^^ includes all the relativistic corrections up to order M^^ i.e. in addition to the 
non-relativistic spin current and convection current terms of order it includes all the terms of order 
Our expression for this operator given below is the internal operator as defined by Eq. ([6]). We also 
assume that the initial momentum Pj is directed along q. (This is the case for the ANB reference frame.) 
Then the current operator includes dependence on q, uj and the magnitude of Pj. In the expression for it 
below all the momentum operators are placed on the right hence rendering the operators in non-symmetric 
forms. Nevertheless the Hermiticity of J(x) is still intact but in momentum space reads as J^(q) = J(— q). 
We use the notation r' = r — R and p' = p — A~^P, P and R being the total momentum operator and 
the non-relativistic center of mass operator. 

The resulting one-body current operator is 

J(i) (q, uo, Pi) = ispin + jp + j, + Aj + (cu/M)j^, (8) 
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with 



J spin 
ip = 

ii -- 

Aj 



t[a X qj 
2M 



M 



G 



E 



2 



2M 



8M2 

,2 



1 - 
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(k' + 2) 



+ 3) 



G 



M 



+ G 



M 



'MP 



(9) 
(10) 

(11) 



+ [Gm-Gz,(1 + 2/€2)] q(p'.q; 
-2zG^[axq][(p')' + «:(p'.q)] 
- Gm)[p' X q] [K{a ■ q) + 2{a ■ p')] |, 



and 



I Ge — 1G 



M 



8M 



^q + iK[cf X q\ + 2i[a x p']) 



(12) 
(13) 



In the above expressions we use the notation 



Gr_m - G^e j^j{Q'^) — ^ + G% m{Q'^) — ^ 



where are the Sachs form factors. We also denote 



(14) 



K = 1 + 2Pi/Aq. 

(Note that 2p + q = 2p' + /tq.) The terms jp and jg together represent the convection current. The latter 
longitudinal component of this current does not enter the net response. However, this component, which 
contains the charge operator, is required when one uses an alternative expression (the Siegert form, see 
Eq.(20) fromRef. R) for electric multipoles based upon the continuity equation. If we chose to use this 
form in the present calculation then the charge operator with inclusion of the standard Darwin-Foldy and 
spin-orbit corrections would be used for calculating electric multipoles. In detail this is given by [23 1 



p(q,u;) 



Gr 



q^ 



Ge — 2(j 



M 



i{a ■ [q X p']) 



(15) 



m 



8M2 J 4M2 

i^'or the two-body current operator we use the customary non-relativistic expressions, of the form listed 
4], Appendix A. The regularization constants entering the two-body current are adjusted to the NN 



interaction we use so that the continuity equation is satisfied approximately, see 



12j. For high q values 



the relative contribution o: 
As explained in Ref. 



the two-body current in the region of quasi-elastic peak is less important. 
J], the current operator is to be used in the form of an expansion over the 



multipole operators Tjj^{q,uj) and T^^{q, 



UJ 



A=el,mag jra 



(16) 



Here e = in the electric case and e = 1 in the magnetic case. The quantities are electric and 
magnetic vector spherical harmonics [l^. We calculate the multipole operators Tj^ in terms of similar 
operators Tj^ related to the vector spherical harmonics of the form 

Y5™(q)= Yl CCYim'{^)e,. (17) 

m'+fi=m 



14| . and I = j ± 1, j . From expressing the expansion of Eq. 



Here are the spherical unit basis vectors 
(fT6|) in terms of the harmonics (fT7j) one obtains the operators 

^i- = 4^ / ^q(Y5^(q) • J(q,^,^.)) • (is) 

These operators are irreducible tensors of rank j. In accordance with the expressions for the harmonics 
^fm^'^^ in terms of the harmonics (fTTl) one has 



f;:r = fj„. (20) 



Expressions for the components of the multipoles fll8p pertaining to the current ([H]) are listed in the 
Appendix. The alternative expression for electric multipoles of the current contains also the multipoles 

Pjmiq) = ^ j dqyjrn{ci)piq,uj) (21) 

of the charge density operator. 



The dynamical part of the calcu 
as for the lab response function in 



ation of the response function R^^^ ((Tj) is performed in the same way 
4|. 



IV. RESULTS AND DISCUSSION 

As mentioned in the introduction we use the AV18 NN potential and the UIX 3NF as nuclear force. 
The calculation is carried out in the ANB frame for eight momentum transfers ^anb: 400, 450, 500, 550, 
600, 650, 700, and 750 MeV/c. We consider electric and magnetic multipole contributions up to a maximal 
total angular momentum J™^^ of the final state such that a convergent result of Rt is obtained for any 
q value. For instance we take JJ^^^ = 19/2 and 37/2 for q = 400 and q = 750 MeV/c, respectively. 
As already poined out before we use the LIT formalism [3, Q] in order to take into account final state 
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interaction. For the LIT parameter o"/ we choose two different values, namely cr/^i = 5 MeV and = 50 
MeV. We combine both results in the following way 

Ltotiaji,ai) = L{an,cTj,,)f{aR) + (^^^ LK, a,,2)(l - /K)) , (22) 

where L denotes the Lorentz transforms of the response, and 

fian) = exp{-{an/aof) (crn > 0) and /(an) = 1 (a« < 0) (23) 

with (To = 100 MeV. This choice has the advantage that one has a relatively large resolution for the Rt 
behavior at lower energies, while for the high-energy behavior a smaller resolution is completely sufficient. 
The inte gra l eq uation that corresponds to the transform Lt^t was solved to extract Rt- The inversion of 



the LIT 



15 



If 



17( 1 has been made as described in 
In Fig. 1 we show i?T (q'anb, i^anb) for the above mentioned eight q values using our full current operator 
(relativistic one-body + isovector MFC consistent with AV18). One sees that the ^anb dependence of Rt 
exhibits a very regular and smooth pattern. This allows us to use a spline interpolation to determine 
-Rt(?anb5 ^anb) for intermediate ^anb values. In this way we are able to obtain results for RT{(liab,^ia.h) 
via the transformation of Fq.(l5]) for 500 MeV/c < giab < 700 MeV/c. 

In the following we investigate three different theoretical aspects: (i) comparison of lab and ANB frame 
calculations, (ii) relativistic contributions to the one-body current operator, and (iii) the MFC contribution. 
We first turn to the comparison of lab and ANB frame results. In Fig. 2 we show -Rr('?iab, "^lab) evaluated 
with the nonrelativistic one-body current for lab and ANB frame calculations. The ANB results show a 
sizable shift of the peak position to lower energies, which grows with increasing q. In detail one has the 
following shifts, 8.7, 16.7, and 29.3 MeV at g = 500, 600, and 700 MeV/c. The size of the shifts is very 
similar to those found for the longitudinal response function Rl in and corresponds to the differences 
of non-relativistic and relativistic kinetic energies of a nucleon with momentum qiab (see discussion of peak 
position in section III). One also finds an increase of the peak heights, namely by 5.6%, 10.3%, and 16.7%. 
The relativistic contribution to the one-body current is illustrated in Fig. 3. It leads to a reduction of the 
peak heights of 6.2%, 8.5%, and 11.3% at g = 500, 600, and 700 MeV/c, while there are no sizable effects 
on the peak position. Finally Fig. 4 shows the MFC contributions. As one might expect they are rather 
small and decrease with increasing q. In detail one has increases of 3.2%, 2.7%, and 2.2% for the three 
considered q values. 

Now we turn to a comparison with experimental data (see Fig. 5). For all the three considered mo- 
mentum transfers one finds an excellent agreement of experimental and theoretical peak positions. For 
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q = 500 MeV/c one also has an excellent agreement of the peak height. At g = 600 and 700 MeV/c 
the theoretical peak height overestimates the data by about 5% and 10%, respectively. We would like to 
mention that a different choice for the nucleon form factor fits should lead to rather small effects only. 
The reason is that at higher momentum transfer Rj- is dominated by the spin current contribution where 
the magnetic nucleon form factors enter which for the various fits are rather similar in the range 500 
MeV/c < q < 700 MeV/c (e.g., compare the dipole fits with those from 18|). In the present work we 
do not consider any A degrees of freedom. As shown in 2|, up to g = 500 MeV/c there are only tiny 
A effects in the quasi-elastic region. Also at higher q one may expect that the quasi-elastic response is 



not affected much by A isobar currents (compare to deuteron electrodisintegration results, see e.g. |l9l]). 
The increasing difference between theory and experiment with growing momentum transfer suggests that 
unincluded relativistic effects (wave function boost, dynamical effects) are increasing in importance. In 
future we will investigate to see if we can get a better understanding of these effects in order to improve 
the comparison with experiment. 

For the comparison with the experimental data of Fig. 5 one has to consider that pion production is not 
taken into account in our calculation. The pion production thresholds are at about 180, 200 and 220 MeV 
at g = 500, 600, and 700 MeV/c, respectively. For q = 500 MeV/c one can nicely see that the theoretical 
Rt starts to underestimate the experimental Rt in the pion threshold region. 

To sum up we can say the following. We have calculated the ^He transverse response function Rriq, uj) 
with a realistic nuclear force (AV18 two-nucleon and UIX three-nucleon potential) in the quasi-elastic 
region at 500 MeV/c < q < 700 MeV/c with full inclusion of final state interaction. The calculation is 
carried out in the ANB frame with a subsequent transformation of Rt to the lab system. Relativistic 
effects to the one-body current operator as well as meson exchange currents are taken into account. The 
relativistic effects reduce the quasi-elastic peak, while the MEC contributions are rather unimportant. 
The use of the ANB frame provides excellent agreement with experimental peak positions. Concerning 
the peak heigts one finds a good agreement of theoretical and experimental results at g = 500 MeV/c, 
while theory overestimates data up to 10% at higher g. 
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APPENDIX A: MULTIPOLES OF THE ONE BODY CURRENT AND CHARGE OPERA- 
TORS 

In the formulae below we use the notation 



ijj = jj{qr')Y,Ur'), = V2^Tl, Uab = ^/ i2a + 1) i2b + 1) . 

The quantity d'^^ below is defined by the relationship —id'^ = p'^, and = {d' ® <9')7/i- Denoting 
_j^Q>{A) _ where p'^ is the last particle internal momentum, one has 

1/2 Q 



A-l 
A 



Here the derivative is taken with respect to a component of the last Jacobi vector defined as ^a-i 
Various operators entering the current ([8]) give the following contributions to the multipoles (fT8l) : 



{Am^-y'l dcie'^^'t (Yj„(q) ■ [a X q]) 



1/2 



J + 1 



{ATTZ^y I dqe^^^'t (Yj+^(q) • [a x q]) = - 
{Am^y [ dcie'^^'z (Yj-^(q) • [a x q]) = - 




{iPj (g) a] 



jm ' 



jm 



{Am^-y I rfqe^'^'-' (Yj„(q) • p') = {^|J, ® 9'),, 
[Am^y I dqe^^^' (Yj±^(q) ■ p') = ± ® d'] 

{Am^-y I dqe^^^' (Yj„(q) . q) = 0, 
(47r^^r^y" rfqe^'^^' (Yj:'(q)-q) = 
(47rz^)-' / t/qe'^^r' (Yj-^(q) ■ q) = 



jm 



J + 1 



2j + 1; 

1/2 



2j + l 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 

(A6) 
(A7) 
(A8) 
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' J dqe^'i^' (Yj^(q) • p') (p' • q 



1 1 2 

J - 1 J J 



1 1 2 

J + 1 J J 



{Attz^)-' I dcie'^^' (Yj:^(q) • p') (p' • q) = 
{Ani^y I dqe^^^' (Yj-^(q) • p') (p' • q) = Sj, 



7=0,2 I A - 2 J A - 1 



7=0,2 I A J A - 1 



{A^V-'Y' j dcie'^^' (Yj^(q) • q) (p' • q) = 0, 
{^nVy' j dqe'^^' {Y^^'ik) ■ q) (p' • q) = VT+l 5, 
{A^Vy' j dqe'<^r' (Yj-^(q) • q) (p' • q) = - v7 S, 

S = (2j + 1)-^ [vjlV'i-l ® d')jm + V7Tl(V',+l ® d')jr. 



X 



^^^^ 

(2j-l)(2j + l) 



{Ani^-y' J dqe'^^^'z (Yj^(q) • [a x q]) (p' • q) 



(2j + l)(2j + 3) 



1/2 r 



(47ri^)-' J dqe^^r'i (Yj;:,^(q) • [a x q]) (p' ■ q) ^ y/j S, 
{iTTi^y J dqe^^^'i (Yj-^(q) ■ [a x q\) (p' • q) = y^^, 



jm 
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= (2j + 1)-' 



jm 



{Am^-') j dqe'^'^'t (Yj^(q) • [p' x q]) (a • q) 

1 7 7 - H ^ /jm 



X ^ (-i)'-^n, 



1 J J 

1 i J - 1 



1 j I 



i j j + i 
1 j I 











jm_ 






1/2 


L 2j + 1 






rj(j + i)" 


1/2 


jm 


L 2j + 1 . 








) jm 



2j + l 



n 



1 J + 2 J + 1 



1 j j + 1 J ^ ^^"^ 
(47ri^) y dqe^'^r'z (Yj+^ (q) ■ [p' x q]) (a • q) = VJ S, 
{Am^y J dqe'^^'i (Yj-^(q) • [p' x q]) (a • q) = ^/j + lS, 

I 3-1 3 



'^3 + '^J I 1 J / 



jm 



+ 



i + 1 



2j + 



1/2 



^ '=j,j+i 1^3 ' 



J dqe^'^^'t (Yj„(q) • [p' x q]) (a • p') 



1/2 



+ 



I ; ■ 

7=0,2 / [ ' i ^ 1 J 

I 1 1 7 



2j + 



7=0,2 



I 3 + i 3 



{{i/jj^i (g) X^)^ (j)^. 



(47ri^)-^ / ciqe^'^-'i (Yj^^lq) ■ [p' x q]) (a ■ p') = ( ^ 



1/2 



(47ri^)-' J dqe'^'^'i (Yj-^(q) • [p' x q\) (a ■ p' 

1 1 7 



J + 1 



1/2 



7=0,2 I \ t 3 3 



2j + l 
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I dqe'^^'z (Yj^(q) ■ [a x p']) 



5Z ^ ^ ^ ^ {i^j®d')^(S)a) 

i=j±i,j I J J ^ 



(A27) 



(Attz^)-' I dqe^^^'^ (Yj^^(q) ■ [a x p']) 



(A28) 



At deriving these formulae the expressions for n ■ Y'„j(n) in terms of spherical harmonics, and for 



nYjm{n.) and [n x Y'„j(n)] in terms of vector spherical harmonics [l^ have been used. We used also the 



relationship 



(n.a)(YWn).b) 



1/2 



((l^_i(n)®a),®6),„- 



/ + 1 



1/2 



{{Yi+i{n)0a)i0b) 



21 + 1 J ' ' 'jrn + 1 

The spin-orbit component of the charge density operator (fTSll in (l2TI) leads to the multipoles 



' j dcie'^''Y,m{(iW " [q x p']) 



2j + l 



1/2 



1 1 1 

J + 1 J I 



(A29) 
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FIG. 1: (color online) Rt{qanb, Eanb) of ^He with relativistic one-body and meson exchange current at various 
q values (internal excitation energy -Eanb = i^anb + Qanb/^t)- 
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FIG. 2: (color online) R-Tisix&h^^x&h) of ^He from ANB (dashed) and lab (dotted) frame calculations with 
relativistic one-body current. 
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FIG. 3: (color online) -RrCo'iab, "^lab) of ^He from ANB frame calculation with relativistic (dash-dotted) and non- 
relativistic (dashed) one-body current. 
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FIG. 4: (color online) i?r('7iab, ^^lab) of ^He from ANB frame calculation with relativistic one-body current with 
(full) and without (dash-dotted) meson exchange current. 
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FIG. 5: (color online) RT{Qia.h,^iab) of ^He from ANB frame calculation with relativistic one-body and meson 



exchange current (full) in comparison to experimental results from 



20(| (squares), ^21|] (diamonds), [22j] (circles). 



19 



